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Abstract. The neutron star's crust and mantel are typical examples of no n- uniform bulk 
systems with spacial localisations. When modelling such systems at low temperatures, as is 
the case in the crust, one has to work with antisymmetrised many-body states to get the 
correct fermion behaviour. Fermionic molecular dynamics, which works with an antisymmetrised 
product of localised wave packets, should be an appropriate choice. Implementing periodic 
boundary conditions into the fermionic molecular dynamics formalism would allow the study 
of the neutron star's crust as a bulk quantum system. Unfortunately, the antisymmetrisation 
is a non-local entanglement which reaches far out of the periodically repeated unit cell. In this 
proceeding, we give a brief overview how periodic boundary conditions and fermionic molecular 
dynamics can be combined without truncating the long-range many-body correlation induced 
by the antisymmetry of the many-body state. 



1. Introduction 

Neutron stars represent one of the densest forms of matter in the observable universe. As an 
"outgrown nucleus", they confine approximately 10^^ baryons inside a radius of roughly 10 km. 
This leads, in the centre, to densities expected to be as large as five to ten times the normal 
nuclear saturation density uq. In the outer layers of a neutron star, the outer and inner crust 
which represents only a tiny fraction of the total mass, the densities are on the average below 
no- Nonetheless, it is associated with various astrophysical phenomena such as the cooling of 
neutron stars or neutron star glitches [Ij. A detailed knowledge of crustal matter will help to 
improve our understanding of these phenomena. 

The neutron star's crust consists of protons, neutrons and electrons that organise them- 
selves as a result of the combined effect of the short-range nuclear attraction and the long-range 
Coulomb repulsion. In the outer crust, i.e. the lower-density regions, nucleons bind to form 
nuclei that are placed in a Coulomb lattice embedded in an electron gas. In the inner crust, 
the neutron-drip density 4 x lO^^g/cm^ [2j) is surpassed giving rise to neutron rich nuclei 
placed on a similar Coulomb lattice, immersed in a gas of neutrons and electrons. In the more 
dense layers of the crust (no/3 to no/2 [3J), the spherical nuclei outgrow the lattice spacing. 
In this crust-core interface, i.e. the mantle of the neutron star, it is believed that the system 
balances on a subtle interplay between the nuclear surface energy and the Coulomb energy of the 
neutron-proton-electron system. This results in a multitude of competing quasi-ground states 



with similar energies, often having quite different matter distributions where the system can 
cool down to. These complex-shaped foam like structures are referred to as "nuclear pastas" 

dims]. 



2. Fermionic molecular dynamics 

The crusts and especially the mantel of a neutron star are excellent examples of bulk fermion 
matter where spatial localisations play an important role. Many computational techniques have 
been used to investigate these. These range from the liquid drop model over Hartree Fock to 
Monte Carlo and molecular dynamics ^41 El [71 [HI O HOl |TT]. With a few exceptions, all of these 
methods study a "single" geometry of the nuclear cluster. This geometry is either explicitly 
implemented or strongly correlated to the shape of the unit cell and the imposed boundary 
conditions. However, the implementational or computational constraints should not distort the 
properties of the macroscopic system. The Monte Carlo and molecular dynamics methods are 
unbiased with regards to the geometry of the nuclear clusters and allow the study of more re- 
alistic systems. As there are many configurations with almost the same energy, the nuclear 
pasta phases are susceptible to small perturbations stemming from external probes or slight 
temperature changes. The molecular dynamics picture is appropriate to study such systems 
and is already successfully used for classical descriptions. From the many existing approaches 
[121 [El [HI 113 [IS], classical molecular dynamics (CMD) provided a way to study the effect of 
large density fluctuations on neutrino opacities and evaluated the breaking strain of the neutron 
star's crust [3 [17]. Quantum molecular dynamics (QMD), on the other hand, identified various 
pasta structures and evaluated one of the first phase diagrams of the different pasta phases [8]. 

Disregarding the classical nature of CMD and QMD, their achieved successes in the study of 
nuclear matter at sub-nuclear densities have given us a first understanding of the crust's bulk 
properties. However, when studying phenomena at length scales smaller then the de Broglie 
wavelength of the system's constituents, quantum effects stemming from the wave character 
of the nucleons predominate. In addition, the statistical many-body correlations related to 
the system's fermionic content will play a fundamental role. While CMD lacks both the wave 
character as well as the statistical description, QMD works with density packets representing 
distinguishable particles. The indistinguishablity of the fermions is mimicked by means of a 
phenomenological Pauli potential which prevents identical fermions to occupy the same phase 
space [18]. Another fundamental quantum feature missing in QMD is the dynamical variation 
of individual spin components and the spreading of the wave packets. 

The shortcomings of QMD for fermion systems are remedied in fermionic molecular dynamics 
(FMD) [151 [Ej- While QMD uses a phenomenological Pauli potential to mimic fermion 
statistics, FMD introduces the fermion many-body states as antisymmetrised products, i.e. 
Slater determinants, of localised single particle states \qp(zp)) : 

|Q(zi, . . . , za)) = A\qi{zi)) • • • \qA{zA))- (1) 

This way, the Pauli exclusion principle and uncertainty principle are implemented from the 
beginning and cannot be violated. With a proper choice of the single-particle states \qp(zp)) 
and their variation parameters z^, concepts as average phase- space position, spin, isospin and 
spreading can be introduced. While the equations of motion of the many-body fermion system 
can be obtained by means of the time-dependent variational principle, structure studies are 
obtained by means of Ritz' variational principle jT9] . 



When studying the observable properties of fermion systems represented by many-body states 
as one has to evaluate the expectation values of various operators. Although the single- 

particle states in the FMD wave function are not orthonormal, the expectation values can easily 
be evaluated by means of a matrix formalism as 



A 1 ^ 

^/ = ^ {qp\Bl\qq)Oqp, mu = 2 {QpQr\13ll\qqqs){OqpOsr " OqrOsp)- (2) 

pq=l pqrs=l 

The wave function's determinant structure is embodied in the matrix o, representing the inverse 
of the overlap matrix n with Upq = {qp\qq)- These two matrices are fundamental for a correct 
fermion description. Therefore, and because their eigenvalues can cover many orders of magni- 
tudes, it is paramount to calculate these matrices as accurately as possible (i.e. analytically). 

FMD was originally devised to describe heave-ion reactions in a time-dependent framework 
[151 [16J. Now, FMD is also successfully used to describe stationary situations related to nuclear 
structure and nuclear reactions by mixing many Slater determinants [191 [20]. However, due to the 
long-range character of the Pauli correlations, evaluating bulk fermion matter becomes a tedious 
task. The dimension of the matrices usually impedes simulations of a large number of particles, 
however desirable these may be for bulk fermion systems such as crustal matter. Applying FMD 
on bulk fermion systems has therefore been a long open question. This has recently been solved 
within the concept of "periodic boundary conditions" by making the spatial positioning of the 
single-particle states periodic [2L,22j. An infinite fermion system is created by replicating a unit 
cell — containing A particles — on an infinite lattice 53. Each cell can be identified by a lattice 
vector R = niai n2a2 -\- n^as with integer nj. Using the translation operator T{R), the FMD 
many-body state can be written as 

\Qoo{zi, . . . , za)) = i (g) f{R) {\qi{zi)) • . . \qA{zA))} • (3) 

As indicated earlier, the overlap matrix n and its inverse o are fundamental for the correct 
fermion description of the system. Although both matrices have infinite dimensions, they 
exhibit a peculiar nested block- Toeplitz structure which can be exploited [21j. As a result 
of the translation invariance of the system, each A x A block of the overlap matrix can be 
identified unambiguously with the lattice vector R that connects the two cells of the bra and ket 
states. The blocks can be evaluated as iipq^R = {qp\T{—R)\qq) . Although n and o are of infinite 
dimension, it is possible to obtain the inverse overlap matrix through the following scheme 

Af{k) = n^e-^^-^, 0{k)=Af{k)-\ OR=-^ [ 0(fc)e^^-^d3fc, (4) 

^BZ JbZ 

where BZ represents the first Brillouin zone of the lattice 53 and is a vector in this volume [2T\ . 

The strength of the proposed formalism reveals itself upon evaluating the operators in 
reciprocal space. Normally, the expectation values of the infinite fermion system would be 
calculated by means of Eqs. ^ resulting in infinite sums over the block structure. These 
sums, however, translate into integrals over the first Brillouin zone which are computationally 
straightforward to evaluate. The expectation value of a one-body operator per unit-cell volume 
can then be computed as 



Vbz 



A 
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Here we assume that the operator Bj commutes with the translation operator T{R) over any 
lattice vector R. The two-body operator has a similar structure. 

3. Results 

When studying properties of bulk matter, it is important that the constraints of the simulation 
do not interfere with the result and that the properties of the simulation should be equivalent 
with those of the macroscopic system. It stands to reason that for non-periodic systems, the 
periodicity should not influence the simulation's results. As long as the correlation volume of 
the interaction does not exceed the simulation volume, the imposed periodicity works fine. How- 
ever, serious problems arise in the presence of long-range correlations such as induced by the 
Pauli exclusion principle. We demonstrate that the proposed technique reproduces two essential 
features intrinsic to free fermions. The single-particle states are Gaussian wave packets of the 
form {x\ab) = exp{ — (cc — 6)^/ (2a)} where the complex vector b represents the mean phase-space 
position and a is a complex parameter connected with the spreading of the wave packet in phase 
space. 

In Ref. [15] it was shown that, for a hundred periodically positioned wave packets, Eqs. ^ 
reproduce the momentum and spatial densities intrinsic to a one-dimensional Fermi system. An 
identical result was achieved in Ref. [21j by using periodic boundary conditions with a unit cell 
containing a single particle. The densities were evaluated using Eqs. In this reference, it was 
shown that with increasing overlap of neighbouring wave packets, the effect of antisymmetrisa- 
tion makes the fermionic behaviour apparent. This effect is also seen in Fig. [T] which compares 
the kinetic energy of the same periodic system with that of an individual Gaussian and that 
of a free one-dimensional Fermi gas as function of the wave packet spreading. Although this 
one-dimensional case hints towards free fermion behaviour, it must be mentioned that, in one 
dimension, the free Fermi gas results and that of a weak periodic system coincide when having 
identical particle densities. In more dimensions, it can be shown that the one-particle-in-a-box 
results evolve towards those of a weak periodic system [21j. This effect is expected because with 
increasing overlap, ergo increasing wave packet spreading, the wave packet seems uniform to the 




Figure 1. The kinetic energy of periodic Gaussian wave packets (solid line) is compared with 
that of a single Gaussian (dashed line) and the free Fermi gas (dash-doted line). The kinetic 
energy per particle is scaled with the mass m of the fermion and the periodic spacing £ 
(box-size) of the Gaussians. The kinetic energy is presented as a function of the scaled variance 
^Mea/i of the Gaussian wave packets. 
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Figure 2. Comparison of the scaled density distributions (left four panels : spatial px^^ /N] right 
four panels: momentum pj^kJ/N) of = 25 Gaussian wave packets under different conditions. 
The wave packets are randomly placed in the unit cell of a square lattice with unit length I 
and have a width parameter given hj a — (0.2^)^. The mean momentum of the individual wave 
packets is set to zero. The spatial densities are normalised and coordinates are expressed in units 
of I. Analogously for the momentum density with coordinates expressed in units of A:^ = \fNi\ ji. 
For each type of density distribution : the upper-two panels depict distinguishable particles, the 
lower-two panels represent the antisymmetrised case, the left and right panels show the effect 
of a calculation without and with periodicity, respectively, as described in the text. For the 
spatial density, the black circles show the positions of the centroids of the wave packets. For the 
momentum density, the dotted circle represents the Fermi "sphere" of a system of free particles 
with the same mean density, — y/^nN ji. Note the difference in colour scaling. 



unit cell and the system behaves as if it consisted of plane waves in a weak periodic system. A 
successful investigation of bulk matter obviously requires that influences of the geometry of the 
chosen boundary conditions are negligible, and thus the wave packets should be smaller then 
the unit cell. 



Fig. [2] depicts the effect of periodicity and antisymmetry on the spatial (left hand side) and 
momentum distributions (right hand side) for 25 single-particle states with zero mean momentum 
randomly distributed in a square unit cell. For distinguishable particles (upper four panels) the 
momentum distribution of the 25 particles is the same as the one given by the individual wave 
packet. When enabling the periodic boundary conditions, only the spatial density grows at 
the boarders due to the tails leaking in from the neighbouring cells. For indistinguishable 
fermions (lower four panels), where the many-body wave function is antisymmetrised, one sees 
that the spatial density (without periodic boundary conditions) in the lower left most panel is 
much smoother. The fermions tend to avoid closeness in coordinate space and at the same time 
enlarge the occupied momentum space as indicated in the corresponding panel for the momentum 
density, thus, inducing Fermi motion. When enabling the periodic boundary conditions, the 
spatial density increases again due to the tails of the neighbouring cells. In momentum space, 
the distribution is even more spread due to the antisymmetrisation which is now also with 
respect to the particles in all surrounding cells, neighbouring and beyond. One obtains an 
almost uniform Fermi distribution where deformations in the surface reflect the clustering of the 
spatial density. Hence, a simulation of clustered densities with the proposed periodic structure 



reproduces intrinsic bulk fermion behaviour, unaffected by the imposed periodic structure which 
leaves no trace in the momentum distribution. This clearly shows that antisymmetry, a long- 
range many-body correlation, should not be ignored or truncated for the bulk description. 

4. Conclusion 

This proceedings gave a short overview how periodic boundary conditions can be implemented 
to study fully antisymmetrised, infinitely extended, inhomogeneous fermion matter by means 
of localised single-particle states. The bulk matter is created by spatially distributing a finite 
number of wave packets in a unit cell which is periodically repeated over space. By treating 
each single-particle state as an individual state, the overlap matrix and its inverse are of infinite 
size but possess a nested Toeplitz structure. This is a feature that can be exploited to reduce 
the computational cost and eliminate the infinite dimensionality of the equations. This allows 
us to study bulk fermion matter by means of a finite number of fermions. The structure of the 
resulting equations to evaluate expectation values, resembles those of a finite fermion system, but 
requires an extra integration over the first Brillouin zone. The later is a translation of the periodic 
boundary conditions. Although the equations only address a finite number of particles in a unit 
cell, they keep track of the Fermi statistics of the infinite system. It is demonstrated that this 
formalism incorporates intrinsic bulk fermion behaviour and reproduces density distributions and 
energies of free fermion systems. Furthermore, the use of localised states makes the technique 
suitable to study inhomogeneous fermion matter, such as the neutron star's crust, by means of 
molecular dynamics or Monte Carlo methods. The next step is to include the forces among the 
fermions to study the expected foam-like pasta structures. 
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